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ABSTRACT 

Fallback disks are expected to form around new-born neutron stars following a supernova explosion. In 
almost all cases, the disk will pass through a propeller stage. If the neutron star is spinning rapidly (initial 
period ~ 10 ms) and has an ordinary magnetic moment (~ 10 30 G cm 3 ), the rotational power transferred to the 
disk by the magnetic field of the neutron star will exceed the Eddington limit by many orders of magnitude, 
and the disk will be rapidly disrupted. Fallback disks can thus survive only around slow-born neutron stars and 
around black holes, assuming the latter do not torque their surrounding disks as strongly as do neutron stars. 
This might explain the apparent rarity of fallback disks around young compact objects. 
Subject headings: accretion disks — pulsars: SN87A — stars: neutron — X-rays: stars 



1. INTRODUCTION 

Following a supernova explosion, some of the ejected mat- 
ter m ay remain bound to the remnant and fall back (Colgate 
HUB. A disk can form from this material if its specific an- 
gular momentum I exceeds the Keplerian value at the surface 
of a newly formed neutron star, Z K = \fGMR ~ 1 .4 x 10 16 cm 2 
s , where M is the mass of the star and R is its radius. Nu- 
meri cal simulations of the pre-su pernova evolution of rotating 
stars ( iHeger. Langer & Wooslevl200Cl) suggest I « 10 16 - 10 17 
cm 2 s" 1 , implying tha t disk formation is possible (see also 
IMineshige et all 19971) . 

iMichel & DesslerJ ( 11981 suggested that fallback disks may 
be an ingredient of radio pulsars, as an alternative model to the 
conventi onal view that pulsars are iso l ated e lectromagnetic 
engines. Lin, Wooslev & Bodenheimer ( 1991) associated the 
formation of planets around puls ars w ith fallback d isks. 
iMever & Mever-H ofmeister ( 1989) and Mever-Hofmeister 
i ll 9921) invoked a fallback disk around the remnant of SN 
1987A to explain the deviation of the observed light curve 
from the theo retical li ght curve for pure radioa ctive decay (see 
also IMineshige. Nomoto & Shigevamall993l). 

More rec ently, fallback disks dChat teriee et afl 120001: 
lAlparll200!l iMarsden et alJl2001at lEksi & AlparJl2003l) have 
been invoked to exp lain anomalous X-ray pulsars (AXPs; 
IWoods & Thom pson 2004]) via an accretion model rather than 
the competing magnetar model (D uncan & T hompson 1992T) . 
This model has rekindled interest in the possible implica- 
tions of fallback disks around radio pulsars. A fallback 
disk assisting magnetic dipole radiation torque can explain 
the ag e discrepancy ([Marsden et alJl200Tbl) and braking in- 
dices (Men ou et alJl2001l) of pulsars, as well as the distri- 
bution of radio pulsars on the P — P diagram ( Alo ar et all 
120011) . iBlackman & Pernal ( 12004 suggested that the jets ob- 
served in t he Crab and Vela pulsars are collimated by fall- 
back disks. Oiao et al. (2003) proposed a fallback disk model 
for periodic timing variations of pulsars. If there are fallback 
disks assisting the magnetic dipole torque, field estimates 
based on pure dipole radiation (fi«6x 10 19 (PP) 1//2 ) would 
be overe stimates. Note that the so-cal led high magneti c field 
pulsars dCionzales & Safi-Harbl rlOOH iGonzalez etall 12004 
McLaughli n et alJ l2003albl) have braking indices less than 
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three, indicating that they are not spinning down purely by 
emitting magnetic dipole radiation. 

A fallback disk (TXu. Wang & Oiaoll2003l Ishi & Xull200l 
has also been invoked for the central compact objects (CCOs; 

iPavlov. Sanwal & Teterjl2004l) in supernova remnants. These 
objects are likely relatives of AXPs. A typical member of this 
small group is the central X-ray point source in Cassiopeia 
A which is just 320 years old and hence likely to possess a 
fallback disk. A fallback dis k model of so ft gamma ray re- 
peaters (SGRs ; fV^oods & Thompson 2004) was invoked by 

lErtan & Alparl (120031) who showed that the X-ray enhance- 
ment following the August 27 giant flare of SGR 1900+14 can 
be modeled as the post-burst relaxation of the inner region of 
a disk which had been pushed back by the flare energy. 

While there is ample theoretical motivation for consider- 
ing fallback disks around young compact objects, the exis- 
tence of these disks has never been confirmed observation- 
ally. There are tight limits on the luminositi es of young sys- 
tems like SN87A i lParkll2004l) and Cas A dChakrabartv et alJ 

2001) where fallbac k disks are espe cially likely. For the lat- 
ter object, Chakrab artv et alJ (120011) showed that the optical 
properties of the source are not similar to those of quies- 
cent low-mass X-ray binaries. Infra red emission is predicted 
llPerna. Hernquist & Naravanl 120001) by accretion models of 
AXPs, and possible infrared counterparts have even been 
identified (see Tam et al. 2004, and references therein) for 
five of the six known AXPs. The observations are, however, 
not consistent with a standard thin disk model. Moreover, 
observations of bursts from AXPs (Gavriil, Kaspi & Woods 

2002) suggest that AXPs are relatives of SGRs for which the 
magnetar model is favored. 

In this letter we address the lack of observational ev- 
idence for fallback disks. We argue that fallback disks 
will be short lived because they would be disrupted by 
the spin-do wn power of the neutron star, either in the pro- 
peller phase dShva rtsmaiJl970atlIlIarionov & Sunvaevl 19751 

iDavies & Pringldll98U iRomanova et alJl2004l) or in the ra- 
dio pulsar phase. In the following section we ca lculate the 
spin-down power (|Pr iedhorskv ll986tlShaham & Tavanill98fit 

iTreves & Boccilll987l) transferred from the magnetosphere of 
a neutron star to the disk in the propeller regime and show 
that for conventional initial periods Pq ~ 10 ms and magnetic 

.moments. ~ 10 30 G cm 3 , as estimated for radio pulsars, the 
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luminosity of the disk will be much larger than the Eddington 
limit for all mass flow rates. In this regime, accretion is highly 
unstable, and the disk is likely to evolve very quickly until the 
inner radius of the disk is outside the light cylinder (see § 2). 
Beyond this point, the neutron star becomes a radio pulsar and 
the star can no longer torque the disk (or vice versa). 

In the context of wind-fed, mass-exchange binaries, 
IShvartsmanl {l970b) argued that a fast-rotating neutron star 
would first become an ejector (radio pulsar), in which the 
disk remains outside the light cylinder. Only after the star 
slows down can inflowing matter penetrate the light cylin- 
der, allowing the propeller stage to commence. In the case 
of fallback disks, a fixed amount of mass is available in the 
disk, rather than a continuous supply of matter as in a mass- 
exchange binary. Therefore, once a rapidly rotating neutron 
star progresses through the initial unstable accretion and pro- 
peller stages to an ejector/radio pulsar stage, it is the end of 
the road for the disk; the system cannot switch back to a pro- 
peller or accretion phase. We elaborate upon these points in 
the following sections. 

2. CONSTRAINTS ON A PROPELLERS 

In order that a fallback disk affect the spin evolution of a 
neutron star, the inner radius of the disk must be inside the 
light-cylinder radius = c/fi where is the stellar angu- 
lar velocity. To within a factor of order unity, the inner ra- 
dius of the disk, j? m , is a pproximated by the Alfven radius 
JPavidson & Ostrikerl 197% : 
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where M is the mass flow rate and /i is the stellar magnetic 
moment. When the disk is inside the light cylinder it sup- 
presses radio pulsar emission. If the inner radius of the disk 
is also inside the co-rotation radius R c = (GM/il 2 ) 1 / 3 , the in- 
flowing matter will accrete and spin the star up. 

Initially, a fallback disk is likely to form in the accretion 
phase. As the M in the disk declines (recall that fallback disks 
are not replenished as are disks in binary systems), the inner 
radius moves out, and when it goes beyond the co-rotation 
radius, the system will enter the propeller stage where the in- 
flowing matter, instead of being accreted, is expelled. Here, 
the neutron star spins down through the interaction of its mag- 
netosphere with the disk. If the torque (angular momentum 
flux) acting on the star by the disk is N, the power trans- 
ferred to the disk by the neutron star will in turn be -UN. 
This power will add to the gravitational power of the inflow- 
ing matter, viz., GMM/R m , enhancing the energy budget of 
the disk. Some of the energy will be used to drive matter from 
the system, provided the velocity of the expelled matter, v out , 
is greater than the escape velocity v esc = \j2GMjR m . The 
outflowing matter will carry away kinetic energy at the rate 
(1 /2)Mv 2 ut . Hence, the net radiative luminosity of the disk in 
the propeller phase is given by 

= (GMM/R m -flN-Mv 2 oat /2 ifv „,>v esc , 
dlsk \GMM/R m -VtN if Vout < Vesc- 

The velocity that the expelled matter can attain depends 
on the coupling of the disk with the magnetosphere; i.e the 
torque, which can be written as N = nR^il^R^M . Here, n 
is the dimensionless torque and depends only on the dimen- 
sionless fastness parameter oj* = Sl/fiitC^m)- Conservation of 



angular momentum requires that 

v om =R m n K (R m )(l-n). (3) 

Using this in equation (0 we find that the luminosity of the 
disk in the propeller regime is 

GMM J 1 -w»n-(l -n) 2 /2 if v out > v eS c, 
1 -w,n 



isk ' 
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For Vout < Vesc the matter expelled by the magnetosphere can 
return to the disk at some radius lar ger than the inner radius. 
As shown by Spruit & Taam ( 1993), the matter will then ac- 
cumulate near the inner boundary and accrete sporadically. 

Many prescriptions have been discussed in the literature for 
propeller torques. We estimate n by assuming that the mag- 
netosphere and the inflowing matter are "particles" colliding 
and transferring angular momentum. If the collisions are elas- 
tic and the moment of inertia of the star is much larger than 
that of the inflowing matter one finds n = 2(1 -cj*), whereas 
completely inelastic collisions would give n = 1 — w*. In both 
cases we assume that the mass of the particle representing the 
magnetosphere is much greater than that of the particle rep- 
resenting the accreting fluid. Let us thus define an "elasticity 
parameter" (3 and write 

«(w*) = (l + /3)(l-^), (5) 

where (3 varies between zero and unity and measures how effi- 
ciently the kinetic energy of the neutron star is converted into 
kinetic energy of expelled matter through the interaction of 
the magnetosphere with the disk. The limiting case of (3 = 
corresponds to a "completely inelastic" interaction in which 
the fraction converted to heat is the maximum, while (3 = 1 
corresponds to an "elastic" interaction in which the rotational 
energy of the neutron star is converted completely to kinetic 
energy of expelled matter without heating the disk. This can 
be seen more clearly when we substitute the torque prescrip- 
tion into equation @: 



Ldisk — 



GMM f l + (l-/? 2 )(io,-l) 2 ifv ou t>v e 
2R m \ 2 - 2( 1 + (3){ 1 - cj*)cj* if v ou , < v e 



(6) 



For f3 = 1, Ldi s k = GMM /2R m which is precisely the lumi- 
nosity of a Keplerian dis k with no torque applie d at the in- 
ner boundary. Note that Chatteriee et al. (2000) employed 
n = 2(1 — op , referr ing to the detailed numerical simulations 
of lDaumerid Jl9 96). This corresponds to the limiting case of 
an "elastic interaction" (J3 = 1) in which no heat is produced 
in the disk by the energy transferred from the neutron star. 
As any fluid process that strips matter from the disk would 
be dissipative, the f3 = 1 limit will never be realized in prac- 
tice. In the accretion regime, the tangential velocity of the 
flow will adjust to the velocity of the magnetosphere, which 
corresponds to (3 = 0. For the torque to be continuous across 
the transition from the accretion to propeller regimes, we need 
to assume j3 = at least for w, » 1. As the dependence of the 
expressions in equation on (3 is weak (so long as (3 is not 
too close to unity), choosing (3 in the range 0-0.5 produces 
very little difference in the results. Note that early propeller 
torques employed in the litera ture were quite ine fficient. In 
the original work of llllarionov & Sunvaevl in"975) the dimen- 
sionless torque was estimate d to be l/u>*. In cont rast, the 
recent detailed simulations of Rom anova et alJ (12004) suggest 
even stronger torques than we employ. 

From equation (j6ji, we can solve for the minimum period 
^mitii below which the luminosity of a propeller disk will ex- 
ceed the Eddington limit, Le = 4irGMc/n es , where we take 
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K es = 0.2 for the electron scattering opacity as appropriate for 
the high metal content of the fallback matter. Taking (3 = 0, 
and defining 9 = L^R m jGMM, we obtain 



2nR 



3/2 



(V26 -l + l) ifv out >v eS c, 
(l/2+y^374)" if ^out ^ ^esc- 



(7) 



Another limiting period P m i n , is obtained by the condition that 
the inner radius of the disk be inside the light cylinder. This 

gives 

2tt 



c 



(8) 



Finally, if the disk is in the accretion phase, then the maxi- 
mum M is given by the Eddington limit: GMM max /R = Lem. 
Figure^shows all these limits in the M-P plane. We see that 
only a limited area of the plane is consistent with the various 
constraints we have discussed. Note that, for a given stellar 
magnetic moment, there is a global minimum period P m ; n be- 
low which no stable disk is permitted for any M. A neutron 
star that is born with a shorter period cannot possibly have a 
stable fallback disk. For fi = 10 G cm 3 , the global minimum 
period is P m { a = 35.2 ms. For a general /i and 3> R c , the 
minimum period can be approximated as 

P mm « 0.036^ s, (9) 

where /X30 = ///lO 30 . 

The propeller regi me with super-critical mass flow rates 
has been studied by Mineshige, Rees & Fabian (1991) who 
argued that the inner radius of the disk will adjust such that 
the luminosity does not exceed the Eddington limit. Our 
work shows that for a neutron star with standard parameters, 
/U = 10 30 G cm 3 , P = 10 ms, there is no way for the disk to 
penetrate the light cylinder and yet have its luminosity be less 
than the Eddington limit. In such a system, we argue that the 
accretion and propeller phases will progress very rapidly, on 
a much shorter time scale than the viscous time of the disk, as 
the accretion and spin-down power cause the inner regions of 
the disk to evaporate. Once the disk moves outside the light 
cylinder, conventional wisdom (Shvartsman 1970b) suggests 
that the disk would be disrupted by the radiative pressure of 
the neutron star whic h now acts like a radio pulsar. Recently, 
Eksi & Alparl (|2005) argued that the transition from the in- 
ner zone dipole field to the radiative zone field can be broad, 
especially if the inclination angle between the rotation and 
magnetic field axes is small, which would imply that the disk 
may be able to find a stable inner boundary outside the light 
cylinder. If this is the case, the pulsar would evaporate the 
disk slowly by the interaction of the pulsar wind with the disk. 
Whether the pulsar will be slowed down to periods where the 
inner radius of the disk can re-penetrate the light cylinder de- 
pends on the complex physics of this interaction. 

3. DISCUSSION 

We have shown in this paper that fallback disks around new- 
born neutron stars with conventional periods and magnetic 
moments would be disrupted by a combination of the accre- 
tion power and the rotational power of the neutron star being 
transferred to the disk. Therefore, even though fallback disks 
may well form, they would not be able to survive for standard 
initial neutron star parameters. 

Our results are general, in the sense that for all propeller 
systems with parameters outside the shaded region in Figure^] 



the disk is either outside the light cylinder or will be pushed 
there rapidly. Thus, not all combinations of the magnetic mo- 
ment and period are allowed. In particular, for a given mag- 
netic moment, there is a minimum allowed neutron star pe- 
riod P m i n given by equation <|9}. If a neutron star is born with 
a shorter period that this, only a very short-lived fallback disk 
is possible. This might explain why there has been no direct 
observational evidence for fallback disks around young neu- 
tron stars even though one expects a fair amount of fallback 
material with angular momentum to accumulate during and 
soon after the supernova explosion. 

The recent interest in the fallback disk model is because 
of the potential application of this model to AXPs. An ad- 
vantage of this model is that the neutron stars in AXPs are 
drawn from the same population as radio pulsars, except for 
one new parameter describing the initial mass of the disk 
(Chatteriee & Hernquist 2000). However, our present re- 
sults indicate that, if fallback disks are to be present around 
AXPs, these neutron stars must have been born with long 
periods and/or small magnetic dipole fields that are not typ- 
ical of_jad^o_rMsars ; _n is interesting that recent stud- 
ies IVranesevic et alJl2004l) based on the Parkes multibeam 
survey imply that 40% of pulsars are "injected" with ini- 
tial periods in the range 0.1 - 0.5 s (as first suggested by 
IVivekanand & Naravanll98ll using a smaller sample of pul- 
sars). For the fallback model to be relevant, AXPs must be 
drawn from thi s sub-population of s low in jected pulsars. 

For SN87A, Oge lman & Alparl ( 120041) assumed that the 
compact star was born spinning rapidly and that its present 
luminosity is due to pure magnetic dipole spin-down. They 
concluded that the observed upper limit on the luminosity of 
the star requires the magnetic field of the putative pulsar to 
be either very small (almost in the millisecond pulsar range) 
or very large (in the magnetar range). In the latter case, the 
pulsar has spun down rapidly after the supernova explosion, 
and so the pres ent luminosity is lo w. The magnetar solution 
branch of Ogelman & Alpar (2004) falls in the region of Fig- 
ure 1 in which a fallback disk is not allowed. Therefore, their 
solution is viable; even if the system had fallback material it 
would have bee n ejected early on. Howev er, the second so- 
lution branch of Ogelman & Alpar (2004) corresponds to a 
small magnetic dipole moment (/i3o ~ 0.01) for the neutron 
star, so that a fallback disk is very likely to survive. The lu- 
minosity from the disk in the propeller phase would then be 
much larger than the mag netic dipole luminosity assumed by 
Ogelman & Alpar (2004). Therefore, the upper limit on the 
magnetic field must be much lower in order to satisfy the ob- 
servational constraint on the luminosity. 

An altern ative explanation for the lack of fallback disks 
is given by iFrver. Colgate & Pintol ( 119991) . who argued that 
accretion flows with heavy elements would be dynamically 
unstable when the recombination of heavy elements com- 
mences, giving rise to a line-driven wind that eventually ex- 
pels all the initially bound material. 

Finally, we note that supernova explosions may sometimes 
lead to black holes. The forma tion of fallback disks a round 
black holes was investigated by Mineshige et al. ( 1997), and 
a fallback disk model for ultraluminous X-ray s ources (ULXs; 
Fabbi ano & Whitdl20 04) was suggested bvlQ d2003h . Rotat- 
ing black holes may transfer energy to their surrounding disks 
(e.g., 12002). just like neutron stars, in which case such 
disks would again tend to reach super-Eddington luminosities 
and be disrupted. However, if the energy transfer is signifi- 
cantly less than in neutron stars, the disk might find it easier 
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to survive. It would be interesting to search for fallback disks 
around young black hole systems. Li (2003) mentions two 
ULXs associated with supernova remnants that could be such 
accreting black holes. 
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FIG. 1 . — First Panel: Allowed regions (shaded) and disallowed regions (not shaded) in the (spin period P)-(mass accretion rate M) plane for an accretion disk 
around a spinning magnetized star with a magnetic moment fj, = 10 30 Gem 3 . The boundary between the two regions is demarcated by various lines corresponding 
to (i) R m = Rl which separates the radio pulsar and propeller phases, (ii) R m = R c which separates the propeller and accretion phases, and (iii) L = Le which 
separates the sub-Eddington and super-Eddington accretion phases. Note that there is a minimum allowed period P mm which is independent of M, and a maximum 
accretion rate which is independent of P. Second Panel: Shows the boundaries between the allowed and disallowed regions in the P-M plane for several choices 
of fj,. The minimum allowed period P m j„ scales with fi as in equation |9}. 



